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ABSTRACT 

This paper considers a given set of fixed order compensators for aircraft pitch control problem. By augment- 
ing compensator variables to the original state equations of the aircraft, a new dynamic model is considered 
to seek a LQ controller. While the fixed order compensators can achieve a set of desired poles in a specified 
region, LQ formulation provides the inherent robustness properties. The time response for ride quality is 
significantly improved with a set of dynamic compensators. 


1. Introduction: 

While designing a feedback control, ride and handling qualities are major performance objectives in aircraft 
control problems. Such objectives are normally achieved by dosed closed loop pole assignment [i]. Preserving 
these closed loop poles ( within the desired regions ) in the presence of perturbations is another requirement 
[2]. LQ problems have inherent stability margins to tolerate unstructured uncertainties. LQ design techniques 
with regional pole constraints have been studied extensively in the literature see [3], and its references ]. 
Similar approach, but with dynamic compensators, have been investigated for automotive applications [4]. 
The compensators given in [5] for aircraft control problem are considered in LQ problem setting. The 
objective of this approach is to improve aircraft ride quality defined in [1]. 

II. Problem Formulation: 

An aircraft model in pitch plane [ with normal acceleration ( nd, pitch rate(g ) and elevator deflection ( 
6e ) as state variables and command input ( ) as control variable], is given by[2]: 
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It is well known that the control law 

u{t)=—R~^b Px{t) + r(t) 


( 2 ) 


minimizes the performance index 

J 

and satisfies the algebraic riccati equation 

AP - PbR-^bP + PA + Q=0 


{j: Qx 
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u Ru]dt 


(3) 

( 4 ) 


Selection of weighting matrices to achieve a controller in equation 2 for exact pole assignment has been 
extensively investigated in reference [3]. Suppose, we choose a set of dynamic compensators given in [5] for 
the control law structure [2] ( see Figure-1 ), then the state equations for the compensators are: 
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Figure 2: Regional Constraints for Aircraft 


II. Simulation Results: 

For F-4 aircraft model at Mach=l.. 5, Altitude =35, 000ft, the system dynamic matrices are given by: 

-0.5162 26.96 178.9 

A = -0.6896 -1.225 -30.38 

0 0 -14 

-M5.6 

6=1 0 
-14 

The matrices A and 6 for the sate vector x{t) ==[a:(f); 2 i(f), ■ 22 (f)]’ are 

A 0 0 

A = I Ai -n 0 

A2 0 - Tai 

-175.6 

0 

6 = I -14 I where, 

0 
0 

Ai = [ 1 0 0 ] 

A 2 = [ -0.6896 (-1,225 + t^) -30.38 ] 

At this flight condition, the short period damping ( Qp ) and frequency ( ) requirements are: 

0.35 < Cp<1.3 


3.29< 'Msp< 11.8 

In complex plane, these constraints impose regional pole constraints shown in Figure 2. 

Table- 1 


Design 



Variables 


Csp 

Q.R 

4.5078 

0.4789 

Q,R 

6.44.58 

0.5316 

Desired 

[3.29, 11.8] [0.35, 1.3] 


( 12 ) 

(13) 
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Figure 1: Control Law Structure [Ref 2] with Filters Fi (s) =^^^and F 2 (s) 


Zy -TyZy^ n, (5) 

i'i - q - -T3Z2 -F T 2 Q (6) 

From the aircraft dynamical equations 1, substituting for q, we have 

Z 2 = T {(132 + T 2 )q + a23^e — T 3 Z 2 (7) 


For the new state vector z(t), 
the state space equations become, 
where, 


A = 


b = 


x{t) = [x(f),^i(0, ’ 2 (f)] 
1(f) = Ax{t) + hu(t) 
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It can be verified that for these dynamic compensators, the system In equation 8 Is completely controllable. 
Thus the control law 


u{t) = —R Px{t) + r(f) 


(9) 


minimizes the performance index 


/‘OO 

J= / {xQx+uRu}dt 


( 10 ) 


and satisfies the algebraic riccati equation 


AP-PbR-^bPAPA + Q = 0 (11) 

With the above formulations, we shall now present the closed loop eigenvalues for various values of the design 
parameters. The design parameters for J are obviously the weighting matrices Q and R. Flowever, note that 
the performance index J is significantly influenced by the other design parameters ri,r 2 ,and t 3 , in addition 
toQ and R.Thenext section presents the simulation results. 
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Figure 3: Time response peits due to step input 


The weighting matrices Q =l 3 and R. =10" as well as the weighting matrices Q =l 3 and R = 10" provide 
the acceptabie dosed ioop poies [ seeTabie 1 ]. 

However, what needs to be observed is the time response plots ( due to step input ) shown in Figure 3. We 
observe that the normai acceieration at the sensor iocation is nonminimal. Moreover, the peak acceierations 
are significantiy reduced with dynamic compensators ( about 50%). 
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